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Experimental results on event shapes obtained within (or related to) the method of Effective Charges are discussed in 
view of measurements of the strong coupling, ce s , the /3-function and non-perturbative contributions to event shapes. 
The data strongly advocate to use of the ECH scheme instead of the conventional MS scheme. 



o 
o 

(N 

c 
3 



(N 
> 

m 

in 
o 
\o 
o 

Oh! 
CD 



X 



1. Introductory Remarks 

1.1. Assorted Theoretical Formulae and Results 

Today the MS-scheme is the de facto standard used for comparisons of QCD predictions as well as experimental 
results. It was introduced for practical reason, as an offspring of dimensional rcgularisation, whereas the method of 
Effective Charges (ECH) and the corresponding rcnormalisation scheme are motivated physically. The ECH method 
has been in some detail discussed in the talk of Maxwell |l|. For completeness here the basic formulae and results 
are assembled in view of the later interpretation of the data. 

ECH has been originally introduced in a clear access, taking the observable itself as perturbative expansion 
parameter (called renormalisation group improved perturbation theory, RGI) is given in 3]. With respect to ECH 
I follow the arguments of 0] which applies the ECH method directly to e + e~data. The results especially apply to 
mean values (or moments) of the event shape distributions which depend on a single energy scale, the centre-of-mass 
energy, Q = y/s, only. 

The general next-to-leading order (NLO) expression for event shape distributions of an observable y reads: 
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Here A and B are the first and second order coefficients of the perturbative expansion. Note that due to the 
normalisation to the total cross-section the integral of Eqn. \l\is normalised to 1. The /i 2 (or alternatively the 
X[± = fj, 2 /Q 2 ) dependence reflects the dependence of the prediction on the renormalisation scheme. In NLO the 
change of renormalisation scheme is fully equivalent to the change of the renormalisation scale. 

Weighting Eqn. ^ with and integrating over the observable then yields the mean value (y) which can be directly 
measured. Normalising (y) to the leading order coefficient A a measurable quantity R is obtained. Such quantities 
are called effective charges. The idea of the ECH method is now to chose (in any order of perturbation theory) a 
scheme such that all higher order coefficients for R vanish. This implies that R coincides with the strong coupling 
extracted in a leading order analysis and explains the name ECH as such observables are directly comparable to 
other observables of this kind as well as to the coupling. In NLO this leads to the following scale: 
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An effective charge R obeys the Renormalisation Group Equation (RGE): 
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As the observables R are directly measurable in dependence of the energy the (3r- functions are measurable quantities. 
The leading coefficients of the /3-functions, /?o = 11 — 2/3 • n./ and p\ — (3\/2(3q are universal (the same as for the 
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coupling), the higher coefficients are observable dependent but free of renormalisation scheme ambiguities. Non- 
perturbative contributions, invisible in perturbation theory, lead to the term with the exponential. K is proportional 
to the factor multiplying a power term in the usual power model expressions. Omitting for simplicity the non- 
perturbative term Eqn. [3] can be integrated yielding: 

2 A R R ^ 1 + pxR J Q \ p(x) x 2 {l + pi x)\ ' W 

where p(R) is defined by the square bracketed term in Eqn.|3] The constant of integration, A#, should have physical 
significance and is related to Aj^g by an exact expression |5|. Thus a s (Mz) can be calculated from A# practically 
without additional uncertainty. 

Besides the derivation of the above fundamental equations the paper Q includes a broad discussion of the extraction 
of A^jg from e + e~data and shows the following: 

• The EC H formalism ... is more general than the adoption of a particular scheme ... It can be derived non- 
perturbatively assuming only that the high energy behavior of the observable is given by the leading order 
perturbative expression. 

• ... higher orders can be split into a predictable contribution (of known "RG predictable" terms oc log fi/Q) and a 
remaining piece containing all the genuinely unknown aspects. ... the choice /i — [Iech removes the predictable 
scatter and provides genuine information on the interesting ... higher-order terms. 

• One can exhibit the relative size of the (truly) uncalculated higher- order corrections for different quantities. 
These corrections are related to how the energy dependence of the quantity differs from its asymptotic dependence 
as Q 2 — ► oo. 

• ... by comparing with data one can test how well the first few perturbative terms represent the observed running. 
Marked discrepancies indicate the importance of higher order terms, or ... non-perturbative contributions. 

These findings imply that the ECH scheme offers exceptional experimental control of higher order contributions by 
comparing different effective charge observables or measuring their energy evolution. 



1.2. Remarks on Experimental Data 

The later discussion is mainly based on precise data of shape distributions obtained at the Z-peak and of their 
mean values in the full range of e + e~ experiments above the b-threshold. Some more recent results [HEilSEl refer to 
the energy range ^20-202GcV and comprise data on events with a radiated hard photon (45GeV < y/s < Mz) [SH- 
AH experimental results correspond to simple e + e _ annihilation and are fully corrected for experimental effects. 

It has been realised that hadron mass effects have a particularly strong influence for the so-called jet mass 
observables. This strong mass dependence has been avoided in 0, Q] by redefining the particle four momenta using 
the so called p-scheme: ((p,E) — > (p, \p\)) or E-scheme: ((p,E) — * (pE,E)). Moreover in 0, extra transverse 
momentum from B-hadron decays was corrected for by using Monte Carlo. This correction is power behaved and 
can have influence on the energy evolution of the observables. A figure displaying the correction is included in 0. 

B-decays also influence event shape distributions. For 2-jet like topologies the distributions are depleted. These 
events are shifted towards 3-jet topologies by the extra transverse momentum from the decay. For observables 
calculated from the whole event (like Bt) the decay from the narrow event side contributes for all values of the 
observable leading to visible effects in the whole range of the observable. In principle, the differences between the 
observables due to B-decays need to be corrected when applying power correction models to shape distributions. 
Plots of such corrections are again included in 
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2. Results on Event Shape Distributions 

2.1. Effects of Changing the Renormalisation Scale 




Analyses of event shape distributions performed by experimental collaborations mainly employed the MS-scheme 
with = 1; information about the explicit use of the ECH scheme [ill 113 , fl3| is sparse. Several experiments 
attempted, however, to determine optimal scales by fitting a s and x^ simultaneously to the data. In view of this, so 
called experimental optimisation, it is instructive to discuss the influence of x M on the prediction Eqn. ^ I n Fig- 
the data on r = 1 — Thrust is compared with the prediction Eqn.nf° r 1^ = 1 and the optimised scale (x M ~ 0.0033). 
Evidently the change of scale leads to a "turn" of the prediction. For x^ — 1 the slope of the data is barely described. 
Other observables show a similar though often less pronounced discrepancy. As a consequence the quality of the fits 
for Xfj=l is, in general, unsatisfactorily bad and the fitted values of a s depend on the interval chosen for the fit. This 
is shown in Fig. e.g. for the Thrust [lif. 

This behavior is basically known since Petra 14] or the early days of Lep and has been observed at small values of 
the event shape observables (in the 2-jet region) where the perturbative 0(or s ) predictions are presumed little reliable. 
It should be noted, however, that due to the normalisation of the event shape distributions 1/cr J da/dy dy = 1. the 
misfit between data and prediction must persist at large y if present at small y. This has indeed been seen Il2j. 

The optimisation of a single scale value for a distribution has been criticised Q as, in principle, the scale is expected 
to be y-dependent (compare Eqn. Q). However, the change of scale expected in the typical fit ranges of the data 
is moderate (compare Fig. EJl . Experimental optimisation therefore presents a fair compromise between theoretical 
prejudice or request and experimental feasibility. 

Results on experimental optimisation are included in 1 1 lL I X 3L fl5L 1 1 6j | . The first reference, which is discussed below, 



comprises 18 observables and is based on high statistics data. Moreover in this analysis the dependence of the 
event shape distributions on the polar angle of the event axis with respect to the beam has been exploited. This 
is advantageous in view of experimental systematics and leads to a large number of statistically independent data 
points. 
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Figure 3: Ratio data over NLLA/0(a?)-fit from Opal 00 

(lines added by the author). 
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Consequences for Matched NLLA 0(a^) Analyses 

The mismatch of the slope of data and 0(a 2 ) prediction has consequences also for the analysis of event shape 
distributions with matched NLLA/0(a;?) predictions. The matching from an experimental point of view represents 
a kind of "averaging" of both predictions where in the 2-jet regime the NLLA part and in the 3-jet regime the 0(a 2 s ) 
part dominates. As in the matching the 0{a 2 s ) prediction is used with a renormalisation scale value the slope 

of the data is, for many observables, imperfectly described. In consequences the same problems as discussed above 
for the pure 0(a 2 ) prediction persist, though diminished by the matching with the NLLA part. In consequence the 
X 2 /Ndf is often in-acceptably bad when applying standard rules. Additional theory errors are introduced to regain 
consistency. As in the 0(a 2 ) case the fitted value of a s (Mz) often shows a marked dependence on the observable. 
To illustrate this in Fig. [31 the ratio of the data and the fitted matched predictions is shown for several observables 
from the concluding Opal analysis on event shape distributions 0,^3- in order to make the discrepancies directly 
visible straight lines are additionally put to the data points. For T, C and P>t a clear increase of the ratio with 
increasing values of the observables is seen corresponding to an increase of a s (> 5%). For Mh/E and j/23 the ratio 
is almost constant, while for Byy the ratio decreases with increasing observable. This more complicated pattern 
correlates with the importance of higher order corrections in the 0(a 2 ) part of the prediction or simply with the 
ratio B/A. The three observables with increasing ratio also show a large ratio B/A > 15, for Mh/E and 2/23 the 
ratio is smaller. For {Bw) it is negative. The mismatch of the slope of the data and the prediction also deteriorates 
the initial \ 2 l-^df- The constructed error bands (dominated by "theory errors") restore consistency of the fits. The 
best stability is obtained for j/23 and justifies the small error obtained for this observable (see Fig. Eland [ItI Hsi ) . 

Still the a s (Mz) values depend on the fit range for some observables and are thus biased. This observable dependent 
biases is even present in the overall Lep combination (see Fig. The (correlated) a s values from T and Bt 

are high, j/23 and B\y are low compared to the average. C is close to the average and seems to form an exception. 

It is likely that the imperfect description of the slope by the matched prescription similarly influences analyses of 
power terms introduced in order to describe the hadronisation. E.g. in the corresponding Delphi analysis '6j the a s 
values for Bt and ph — Mj^/s (E-scheme) are markedly smaller compared to the results from Byy, T and p s . 
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2.2. Results from Distributions on Experimentally Optimised Scales 




«; (M z 2)0.1 0.12 0.14 0.1 0.12 0.14 



Figure 4: A\ 2 and a a (Mz) obtained from Oioj) fits of 18 Figure 5: a a (Mz) results obtained from 0(a 2 a ) fits of 18 event 
event shape distributions as function of |l2[. shape distributions for x M =l and a; M =a;° pt . 



In Fig. 0] the A\ 2 and ce s (Mz) as function of as obtained from several event shape distributions is shown 
[ill Il3 |. The fit quality for these fits is satisfactory when including experimental and hadronisation errors only. 
Marked minima in A\ 2 are observed scattered over a wide range of x^ . The clear minima imply that the optimised 
scale values are statistically well determined. The a s {M%) results corresponding to these scales show a much smaller 
scatter as for fixed x^—l. This is more clearly shown in Fig. [5] The spread reduces from 9.4% for x^ = l to 2.2% 
with optimised scales, x° pt . 

It is particularly interesting to study the correlation of the fitted i° pt with the values expected from ECH, x^ CH , 
shown in Fig. EI The x° pt -errors are from the fit, the error bars for x^ CH indicate the range of scales expected within 
the fit interval of the data. A significant correlation p = 0.75 ± 0.11 is observed. The wide spread of experimentally 
optimised scales which has been often criticised is in fact slightly smaller than the range expected from ECH. 

Regarding the distribution of individual observables (indicated by the letters in the plot) it is seen that observables 
calculated from the whole event (e.g. the Thrust) tend to populate the region of small a;^ CH , observables sensitive 
mainly to the wide side of the events (e.g. ph) populate the center while observables from the difference of the wide 
and narrow side show large scales. This ordering corresponds (via Eqn.|2J to large (for small arfj ) and small second 
order MS contributions (for big cc^ CH ). Overall the observed large range of x° pt can be considered as understood in 
ECH theory 

It may be worth noting that the RMS spread 1 of x° pt with respect to xJf CH is only 2.5 (i.e. similar to the range 
used for naive MS theory error estimates). This may be taken as a justification of an ECH theory error estimate by 
a corresponding scale change. However, the corresponding a s uncertainty will be far smaller than in the MS case as 
the optimal scales often correspond to the minima of the a s (Mz)(x f _ l ) curves (compare Fig. 0}. 



1 Estimated from the data of Fig. ©as lo CT(p(lg(a; M /^JT*))). 
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lg (x^ECH) 

Figure 6: Correlation of the experimentally optimised scale a;° pt values with the ECH expectation x^ CH . The errors of x^ pt 
are fit errors, for x^ CH they indicate the expected change in the fit range. 

2.3. Results from Event Shape Mean Values 

Some relevant data on mean values of shape distributions as function of the centre-of-mass energy Q is shown in 
Fig. [3 The data is compared to fits of Eqn. 0] shown as full lines. The dashed lines represent the MS-prediction 
for the same a s (Mz). It turned out that the data can be well described by the ECH/RGI fits with a single value 
of a s (Mz) (~ 0.120) and negligible additional non-perturbative terms. In Fig. [HI the a s (Mz) results obtained from 
this data with MS prediction and Monte Carlo or power hadronisation corrections and ECH/RGI theory with and 
without power terms are compared. Evidently the spread of the results is far smaller for the ECH/RGI results. The 
inclusion of power terms leaves the spread almost unchanged but lowers the average value of a s (Mz) to 0.118. This 
result implies that the differences due to non-trivial uncalculated higher order terms or non-perturbative terms is 
small, of (9(2%). Note that higher order corrections differing by about a factor 2 (in MS) are expected 2 for event 
shape observables calculated from the full event or the wide hemisphere of an event only. This is also reflected e.g. 
in the perturbatively calculated ratio of the power terms of Thrust and ph- The bigger differences seen in the MS 
compared to the ECH results then should be due to "RG predictable" log/i/Q terms. 

The better agreement of the ECH result is also illustrated by Fig. fTO| where the measured size of the power terms 
for some observables as parameterised by «o is compared to the ECH/RGI expectation which has been calculated 
by setting the MS+power terms expression equal to the ECH prediction. It is evident that the RGI/ECH prediction 
describes the data better than a universal value of ag presumed by power correction models. 

Besides the comparison of the a s results from different observables the energy dependence of these observables i.e. 



2 Compare e.g. the B/A ratios in Fig. 14 of 'Qj. 
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Figure 7: RGI fits to several event shape means as function of the centre-of-mass energy (full line) . The dashed line represents 
the MS-expectation (with the same a s (Mz)). 
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Figure 8: q s (Mz) results obtained from several event shape means or integrals over the JCEF or EEC using MS theory plus 
Monte Carlo or power hadronisation corrections and RGI theory with and without non-perturbative terms. 



T002 



7 



FRIF Workshop on First Principles Non-Perturbative QCD of Hadron Jets 



T 
M H 
C 

B w 
2/23 



All 



-+•+- 



I hH- 



J I I I I I I I I I I L 



+•+- 



J I I I I I I I I I L 



0.110 



0.115 



0.120 



0.125 



0.130 

a S (M z ) 



Figure 9: Combined Lep a s (Mz) results for different observables [!£ 
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Figure 10: Size of the non-perturbative terms Qo 
as determined from power model fits compared to 
ECH expectation (grey bands). 



the observed /3r- functions can be compared to the NLO expectation. As the NLO terms present only a small 0(4%) 
correction pi has been set to the QCD expectation. Then the /3 values obtained from the different observables can 
be directly compared. In order to allow full control of the systematic uncertainties this comparison was made for 
seven observables determined from Delphi data only. The measured Po values agree among each other within part 
of the correlated statistical uncertainty and with the QCD expectation (3$ = 7.66 (see Tab. 1). The possible influence 
of power terms (~ 2% at the Z) on the /3-function is small as the energy dependencies of the power terms and the 
/3-function are similar. 

Especially in view of a measurement of the /3-function the analysis has been repeated including reliable low energy 
data on the Thrust. The resulting fit is shown in Fig. ^2 and corresponds to 

/3 = 7.86 ± 0.32 . (5) 

This measurement of the /3-function is the most precise presented so far and allows to strongly constrain the QCD 
gauge group to SU(3) in combination with measurements of the multiplicity ratio of gluon and quark jets [2l| and four 



Table I: Results on f3o obtained from RGI/ECH fits to Delphi data |6J. The first error is statistical, the second is systematic 
and the third is due to the B-hadron decay correction. The uncertainty due to possible non-pert, contributions is small. 



Observable 


A) 


X 2 /N df 


(1 - Thrust) 


7.7 


± 


1.1 ±0.2 


± 


0.1 


9.2/13 


(C-parameter) 


7.8 


± 


1.0 ±0.3 


± 


0.1 


7.2/13 


(p h ) E-def. 


7.5 


± 


1.5 ±0.2 


± 


0.0 


8.8/13 


(p s ) E-def. 


7.5 


± 


1.1 ±0.2 


± 


0.0 


7.1/13 


{B w ) 


7.7 


± 


1.4 ±0.1 


± 


0.1 


6.3/13 


(S T ) 


7.7 


± 


0.9 ±0.1 


± 


0.1 


5.9/13 


(Major) 


8.0 


± 


1.1 ±0.1 


± 


0.1 


9.2/13 
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Figure 11: ECH/RGI fit to data on (1 - Thrust}. The data 
is corrected for the small influence of B-hadron decays. The 
full line represents the QCD expectation, the dashed-dotted 
line the expectation for QCD plus light Gluinos. 
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Figure 12: Constraints on Ca/Cf and Tr/Cf from mea- 
surements of the /3-function, the multiplicity ratio in gluon 
to quark jets and four jet angular distributions. 



jet angular distributions |l9u22l . r2.'j| (see Fig. 112(1 . This agreement in turn lends further support to the measurements 
using the ECH scheme. 

As a final check it has been tried to infer the possible size of pi from the overall Thrust data by fitting the 
corresponding expression with /3q and p\ fixed to the QCD expectation. The resulting value of p2 was found to be 
small and consistent with 0. 

The above quoted results imply that the running of e + e _ data already in the energy range above the b-threshold is 
well described by the leading two coefficients of the QCD /3-function expansion. Neither important NNLO contribu- 
tions nor important non-perturbative terms are observed from the running of the observables or from the comparison 
of a s obtained from different observables. In turn this gives confidence to the ECH/RGI measurements of a s (Mz) 
from event shape means. 



3. Summary 

The extraction of a s (Mz) from seven event shape distribution mean values (corrected for mass effects) using the 
ECH/RGI formalism shows a far better consistency between the individual results (spread ~ 2%) compared to the 
ones obtained using the conventional MS 0{a 2 s ) analysis combined with Monte Carlo or power model hadronisation 
corrections (spread ~ 7%). Moreover the energy evolution of these observables, the /3f?,-functions are well described 
by the two-loop QCD expression and a universal value of a s (Mz). Possible non-perturbative contributions turn out 
to be small (~ 2% at the Z). A fit to the precise data on Thrust in the energy range 15-205GeV precisely confirms 
the two loop QCD expectation. 

These experimental result support the theoretical prejudice that the ECH scheme is more general than a partic- 
ular renormalisation scheme as it can be derived non-perturbatively. Moreover large logarithmic terms oc log/i/Q 
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are avoided in ECH. This fact and the possibility to directly measure the /^-functions and thereby to judge the 
importance of genuine higher order and non-perturbative terms single out the ECH scheme for experimental tests. 

The inciting results obtained for mean values are corroborated by studies of event shape distributions with ex- 
perimentally optimised scales. These scales turn out to be similar to the ECH expectation. Also here the spread of 
the a s (Mz) results is strongly reduced although additional logarithmic dependencies on the observables are to be 
expected. 

In studies employing MS theory and the so-called physical scale 3^=1 the slope of the distribution is often badly 
described. This leads to a dependence of the extracted a s (Mz) value on the fit range, thus to a biased result. Due 
to the normalisation of the event shape distributions to one such a dependence must be present, though to reduced 
extent, in the state-of-the-art matched 0(a^)/NLLA analyses and most likely causes part of the spread of the a s (Mz) 
results extracted using this prescription. Despite the technical problems encountered [l| studies to match ECH and 
NLLA calculations should therefore proceed. 

In the concluding talk of the workshop |3| the question was raised why both ECH and MS /power models lead to 
a reasonable description of the energy evolution of mean values. Taking the ECH results for granted this question 
can be answered: The MS perturbative part of the power ansatz will fulfill the RGE. To leading order this is also 
true for the power contribution as its energy dependence is similar to that of the coupling. As the size of the power 
term is fit to the data the discrepancy to the ECH expectation must be small. Still the predictability of the size of 
the power terms (see Fig. 110(1 using a single value of a s (Mz) as input gives preference to the ECH description for 
mean values. 

It is likely that a large part of the annoying and unsatisfactorily large spread of the a s (Mz) results obtained from 
different event shape observables using standard analyses is due to perturbing logarithmic terms induced by the 
choice of the so-called physical scale and the MS renormalisation scheme. Given the above quoted positive ECH 
results experimentalists should be encouraged to also analyse their data using the ECH scheme where possible and 
not regard this alternative to the MS convention as a heretical digression. 
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